Quantum phase transitions in multileg spin ladders with ring exchange 
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Four-spin exchange interaction has been raising intriguing questions regarding the exotic phase transitions 
it induces in two-dimensional quantum spin systems. In this context, we investigate the effects of a cyclic 
four-spin exchange in the quasi-lD limit by considering a general A^-leg spin ladder. We show by means of a 
low-energy approach that, depending on its sign, this ring exchange interaction can engender either a staggered 
or a uniform dimerization from the conventional phases of spin ladders. The resulting quantum phase transition 
is found to be described by the SU(2)Ar conformal field theory. This result, as well as the fractional value of 
the central charge at the transition, is further confirmed by a large-scale numerical study performed by means of 
Exact Diagonalization and Density Matrix Renormalization Group approaches for A^ < 4. 



PACS numbers; 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 



The destruction of Neel magnetic ordering at zero temper- 
ature by quantum fluctuations has attracted much interest in 
the past two decades.' On top of the possible formation of 
spin-liquid phases, strong quantum fluctuations in spin— 1/2 
systems might lead to unconventional emerging quantum crit- 
icality. In this respect, it has been proposed that there is, 
in two dimensions, a direct continuous quantum phase tran- 
sition between the Neel antiferromagnetic phase and a va- 
lence bond solid (VBS) phase which breaks spontaneously the 
translation symmetry.- Such transition is beyond the conven- 
tional Ginzburg-Landau- Wilson paradigm since direct transi- 
tions between phases whom order parameters enjoy different 
symmetries are generically first-order The exotic quantum 
critical point is described in terms of fractionalized quanti- 
ties that are confined in the Neel and VBS phases and be- 
come deconfined at the transition. The resulting exotic quan- 
tum phase point has thus been dubbed deconfined quantum 
criticality.^ The existence of this new class of quantum crit- 
ical points has been extensively investigated numerically in 
several quantum spin models in two dimensions. "*"" Most 
of these models contain four-spin interactions on top of the 
usual spin-exchange interaction. For instance, the most stud- 
ied model is the so-called J—Q model with a special four-spin 
exchange interaction which is free from the sign problem. "* 
Quantum Monte Carlo computations in this case support the 
existence of a deconfined quantum critical point between Neel 
and VBS phases.^ However, the situation is still controversial 
since other studies suggest a weak first-order phase transition 
which occurs for larger system sizes. ^ 

To get more insight into this problem, we investigate, in this 
paper, the nature of quantum phase transitions driven by four- 
spin exchange interactions in the quasi- ID limit where power- 
ful techniques can be used. In this respect, we consider a iV- 
leg spin ladder with a ring-exchange interaction with Hamil- 
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FIG. 1. (Color online) A-leg spin ladder with leg coupling J^^, rung 
coupling J± and a ring-exchange interaction K. 
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where Sa,n is a spin— 1/2 operator at the n*'' site of the lad- 
der with leg index a = I, . . . ,N. The parameters J|| > 
and J± are respectively the intrachain and the interchain ex- 
change couplings. The ring exchange operator P4 is defined 
on each plaquette [Sa,„, Sa,„+i, Sa+i,„+i, Sa+i,„] between 
two consecutive chains and it cyclically permutes the states of 
the four spins on the plaquette (see Fig. 1 and Eq. (AI)). 

This interaction appears in higher-order corrections in 
the strong-coupling expansion of the half-filled Hubbard 
model'^''' and is known to play a significant role in '^He ad- 
sorbed on graphite.'"*'^ The relevance of the four-spin cyclic 



exchange has also been reported '^^"* in the frame of inelastic 
neutron-scattering experiments for cuprates such as La2Cu04, 
La6Ca8Cu2404i, La4CaioCu2404i, CaCu203, and SrCu203. 

The model (1) in the N — 2 case has been studied exten- 
sively over the years by means of different analytical and nu- 
merical approaches. ^^"^^ The zero-temperature phase diagram 
is rich and several exotic phases have been identified such as 
a scalar chirality phase which spontaneously breaks the time- 
reversal symmetry.^* When J± > 0, the ring exchange desta- 
bilizes the well-known rung-singlet (RS) phase of the stan- 
dard two-leg spin ladder and a staggered dimerization (SD) 
phase emerges. The resulting quantum phase transition was 
predicted to be continuous and to belong to the SU(2)2 uni- 
versality class with central charge c — 3/2.-''''*" The latter 
result was confirmed numerically by means of Exact Diago- 
nalization (ED) and Density-Matrix Renormalization Group 
(DMRG) techniques."''-^^ 

In the general N case, we are not aware, to the best of 
our knowledge, of any investigation of the phase diagram of 
model (1). In the absence of the ring-exchange, it is well- 
known that the physics of the A^-leg spin ladder strongly de- 
pends on the parity of N'^^''^^ . When N is odd, a gapless phase 
with central charge c = 1 is stabilized for all signs of J± 
while a RS phase and a non-degenerate gapful phase, which is 
equivalent to the Haldane phase of the Heisenberg chain with 
spin S = N/2, appear in the N even case respectively for 
J± > and J± < 0."*'"*^ As it will be shown in this work, the 
main effect of a weak four-spin exchange interaction is to in- 
troduce dimerized phases (staggered or uniform dimerization 
depending on the sign of K) on top of these phases. A central 
question is then the nature of the quantum phase transition 
driven by the ring-exchange interaction for N > 2. Can an 
exotic transition, like in the N = 2 case, be stabiUzed? The 
N odd case is intriguing since, on general grounds, one ex- 
pects a Berezinski-Kosterlitz-Thouless (BKT) quantum phase 
transition between a standard c = 1 gapless phase to a dimer- 
ized one as in the Ji — J2 spin— 1/2 Heisenberg chain. "^^ In 
this respect, is it possible to stabilize a more exotic transition 
with extended quantum criticality stemming from the special 
character of the four-spin exchange interaction? Finally, in re- 
lation with the deconfined quantum criticality paradigm, one 
may pose the question of the nature of the quantum phase tran- 
sition when A^ — > 00 to reach the two-dimensional regime. 

In this paper, we investigate these questions by means of a 
low-energy approach when the chains are weakly-coupled and 
ED and DMRG computations for moderate couplings. We re- 
veal the emergence of an exotic quantum phase transition in 
the SU(2)Ar universality class with fractional central charge 
c = 3N/{N + 2) for iV < 4. On the other hand, we ex- 
pect, from the low-energy approach, a weak first-order transi- 
tion for higher N. The SU(2)jv universaUty class has already 
been found in multicritical points of Heisenberg spin chains 
with polynomial interactions,'*'* and more recently in a spin 
chain with three-site interaction.'*^ Similar quantum critical- 
ity appears also in the context of quantum Bose gases with 
long-range interactions.'*^''*^ The SU(2)jv criticality has also 
been shown to describe Haldane-Shastry type spin chains with 
longer range interactions.'*'^''*^ 



The rest of the paper is organized as follows. The low- 
energy approach of model (1) in the general TV case is pre- 
sented in Sec. II. The emergence of the SU(2)7v universal- 
ity class is obtained through two different methods: a semi- 
classical analysis and a more quantum approach based on the 
bosonization technique.'**'*^ In Sec. Ill, we investigate numer- 
ically the quantum phase transition for N = 2, 3, 4 by means 
of ED and DMRG^" computations. Finally, our concluding 
remarks are given in Sec. IV and some technical details on the 
continuum limit are presented in the Appendix. 



II. LOW-ENERGY APPROACH 

In this section, we describe the low-energy approach of 
model (1) in the weak-coupling limit for general N. We 
assume that the in-chain exchange interaction J|| is antifer- 
romagnetic ,/|| > and the chains are weakly coupled: 

|J±,/^|« Jii. 



A. Continuum limit 

In the absence of the interchain coupling (J± = 0) and 
ring-exchange interaction {K = 0), the model boils down 
to N decoupled spin— 1/2 SU(2) Heisenberg chains. As it 
is well-known, the latter model displays a quantum critical 
behavior in the SU(2)i universality class with central charge 
c = 1 or equivalently with one bosonic gapless mode.'*' Using 
the non-Abelian bosonization approach, one can have a con- 
tinuum description of the lattice spin— 1/2 operators in terms 
of their uniform and staggered parts:'** 



Oq ^a,n — Ja ' \ Ij ^a: 



(2) 



where x = nao (uq being the lattice spacing). The uni- 
form part, Ja = JaL + JaR, cxprcsscs in terms of the left 
and right SU(2)i currents JaL.R which generate the underly- 
ing conformal symmetry of the spin— 1/2 SU(2) Heisenberg 
chains. '*''*'*'^' In Eq. (2), the staggered magnetization density 
of the a*'* chain takes the form: 
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Trao 



V2 



Tr (gaa) , 



(3) 



where A is a non-universal real constant, whose value has 
been obtained in Refs. 52-54, and a stands for the Pauli ma- 
trices. In Eq. (3), ga {a = 1, . . . , A^) is the SU(2)i Wess- 
Zumino-Novikov-Witten (WZNW) primary field. This field 
transforms into the spin— 1/2 representation of SU(2) and has 
scaling dimension A = 1/2.'**'^^'^^ It is useful to recall that 
Tr (,9af?) is anti-Hermitian since ga is an SU(2) matrix which 
explains the i factor in Eq. (3). In this continuum description, 
the non-interacting Heisenberg Hamiltonian of each chain ex- 



presses in terms of the currents: 



Hq 



aol 



Y. fdx {jIl + jIb) 

N 
2_^ j dx 3aL ■ 3aR, 



(4) 



where v = 7rJj|ao/2 is the spin velocity and 7 > 0. The 
current-current interaction is thus a marginally irrelevant con- 
tribution which gives rise to the well-known logarithmic cor- 
rections of the SU(2)i quantum criticality of the antiferro- 
magnetic spin— 1/2 Heisenberg chain. ^^'^^ The last impor- 
tant operator in the continuum limit of model (1) is the spin- 
dimerization field of the a*'' chain: 



(-l)"Sa ,„ • S 



A 



a,n+l 



Trao 



V2 



Tr(ga), (5) 



which has the same scaling dimension A = 1/2 as the stag- 
gered magnetization field. It is useful to observe, from this 
representation, that the one-step translation symmetry Ta^ 
corresponds to a Z2 symmetry for the WZNW field: ga — > 

-9a- 

With these results in hand, we can derive the continuum 
limit of the A^-leg spin ladder with ring-exchange interaction. 
The technical details are described in the Appendix for a A^- 
leg spin ladder with general four-spin exchange interactions. 
The interacting Hamiltonian separates into two different parts: 



T~Lint — T~L/\=i + T-Lc 



(6) 



where 'Ha=i is the leading contribution with scaling dimen- 
sion A = 1, which involves the staggered magnetization and 
dimerization fields: 



'Ha=1 = "Y dx {XlUa ■ Ua+l + M eaEa+l) , 



(7) 



with Ai = oo Jl and A2 = SGaoA'/Tr^ in the case of a ring- 
exchange interaction (see the Appendix). The second part of 
Eq. (6) is marginal and expresses in terms of the currents: 



N 



T-icc = Xs^ dx 3 ah ■ JaR 
a=l-' 
N-1 „ 

+ A4 ^ j dx {JaL ■ Ja+lfl + Ja_R ' Ja+li) , (8) 

where we have neglected all chiral interactions and thus ve- 
locity renormalization terms. The identification of the cou- 
pling constants is detailed in the Appendix and we find: 
A3 = ao [-7 + 2i^ (1 - A^) + iX^K (3A2 - l) /tt^] and 
A4 == ao (J_L +4:K - 4A'(1 + AX'^) / n'^) . It is important to 
observe that a A^-leg spin ladder with more general four-spin 
exchange interaction will have the same continuum descrip- 
tion in the weak-coupling limit. The actual form of the four- 
spin exchange interaction is encoded in the expression of the 
coupling constants Ai.2,3.4 (see the Appendix). This means 



that our conclusion on the quantum phase transition for the 
spin ladder (1) still holds for a more general interaction. 

The strategy is to (i) analyze the effects of T-La=i, which 
is a strongly relevant perturbation, (ii) determine the different 
phases of the model and (iii) discuss the nature of the quantum 
phase transition. The next step is then to investigate the stabil- 
ity of these results with respect to the marginal current-current 
interactions (8). 



B. Semiclassical approach 

We will first consider a naive semiclassical approach of 
model (7) which gives some hints about the possible phases 
and transitions of the problem. A more quantum approach 
will be described in the next subsection to further justify the 
results obtained within the semiclassical approach. 

Let us denote W{g) the action of the SU(2)i WZNW model 
which describes the SU(2)i criticality of the spin— 1/2 SU(2) 
Heisenberg chain. This action reads as follows :^^'^^'^^ 

^(5) = ^ / ^'^ Tr(a^5ta^^) + r(g), 



r(.9) 



127r 



d^y e^^'<Tr{g^do.9g^d0gg^d^g), (9) 



where T{g) is the famous WZNW topological term. The ef- 
fective action of model (7) can be obtained using the defini- 
tions (3,5). In this respect, it is convenient to write the product 
of staggered magnetization operators as follows: 



A2 



n-a • "a+1 



27r2a^ 



Tr(,9a'?)Tr(,9a+i'?) 



27r2a2 



Tr(.ga'?)Tr(.gT^iCT 



(10) 



where we used the fact that Tr (^a+if?) is anti-Hermitian for 
an SU(2) matrix. Using the completeness relation cr^^cr* ^ ~ 
2((5q,5(5/j^ — SapS^s/'^), the action of model (7) then takes the 
form: 



JV 



S = Y,W{ga)+ / d^x{Vi+V2) 



a=l 



(11) 



with 



A A^ 

^1 = -^ Y. Tr(.9a.gl+i) 



a=l 



N-l 



V2 - (-Ai + As)-^^ Y. Tr(5a)Tr(g,+i). (12) 

The two contributions of Eq. (12) are of different nature. In 
particular, Vi is invariant under an SU(2)l x 811(2)^; symme- 
try- 9a -^ UgaV, [/ and F being independent SU(2) matrices. 
In contrast, V2 is only SU(2) invariant: ga — > UgaU^- 

We will now apply a semiclassical approach to investigate 
the nature of the zero-temperature phases of model (11). Vi is 
a strongly relevant perturbation with scaling dimension A = 1 



and a gap will open for some degrees of freedom. The main 
hypothesis of the semiclassical approach is to assume that Vi 
operator gives the largest gap Ai of the problem. This hy- 
pothesis is expected to be valid in some parts of the phase 
diagram of the lattice model (1). The numerical simulations 
of Sec. Ill will shed light on the correctness of this semiclas- 
sical approach. Once this assumption has been made, the next 
step is then to write down an effective action which captures 
the low-energy properties of the model when _E <C Ai. The 
nature of the effective action turns out to depend on the sign 
of Ai = ao Jl. Next we discuss the two cases Jl > and 
J± <0. 



1. Jx > 

Since Ai > 0, the configuration ga+i = { — ^)°'9a (a = 
1 , . . . , A^ — 1) minimizes the potential Vi (12) over SU(2) ma- 
trices. Averaging out the .92,. ..,w fields, we obtain an effective 
action on the gi field: 



Seff=NWigi)-N{-Xi+X2) 



A2 



27r2a2 



d'x (Tr(gi))^ 



(13) 

where we used W{—g) = W{g). The latter model is the 
SU(2)Ar WZNW model perturbed by (Tr(gi))^ = Tr($(i)), 
^'^^^ being the spin— 1 WZNW primary field with scaling di- 
mension 4:/{N + 2). This model corresponds to the effective 
action for the spin— A^/2 Heisenberg model obtained by Af- 
fleck and Haldane."*"* The IR property of model (13) depends 
on the sign of Ai — A2 and the parity of N. If Ai — A2 > 0, 
i.e. for a weak four-spin exchange interaction, the minimum 
condition on the potential of the action (13) corresponds to an 
SU(2) matrix with the constraint: Tr{gi) = 0. Since for an 
SU(2) matrix, we have the decomposition: gi = uqI + ia ■ n 
withriQ+n^ = 1, the constraint Tr(gi) = Ogives: gi = ia-n, 
n being an unit vector. Plugging this expression into Eq. (9), 
one obtains the non-linear sigma model with a topological 
term 6 = A^tt.'*'* When N is even, we have a massive non- 
degenerate phase which corresponds to the well-known RS 
phase of the A^-leg spin ladder for J±^ > 0. When N is odd, 
the non-linear sigma model with topological term 9 ~ tt is 
known to be gapless in the SU(2)i universality class. ^°'^' One 
recovers the standard gapless phase of the spin ladder with an 
odd number of legs. 

When Ai — A2 < 0, i.e. K > tt^J±/36, the minimum con- 
dition on the potential of the action (13) is gi = ±/ for all N. 
One enters a dimerized phase (Tr(ga)) 7^ which is two-fold 
degenerate and breaks spontaneously the one-step translation 
symmetry Ta„ (ga — > —go)- The dimerization pattern in a 
given chain is out-of-phase with the one in the neighboring 
chains: (Tr(gi)) = -(1^(52)) = ... = (-l)^-i(Tr(gjv)> 
so that ((-l)"Si.„ • Si,„+i) = -((-1)"S2,„ • S2,„+i) = 
• • ■ = (-l)'^"'((-l)"SAr,„ • ^N.n+i) + 0. This phase cor- 
responds to the SD phase which exists for all N . 

The quantum phase transition between the RS phase, or the 
gapless phase when N is odd, and the SD phase does not de- 
pend on N . It is located at Ai = A2, i.e. Jl = 'i^Kj-n^ . 



The action which controls the transition is 5,;// = NW{g\) 
and takes then the form of the SU(2)Ar WZNW model. We 
thus expect a quantum phase transition which belongs to 
the SU(2)Ar universality class with fractional central charge 
c = 3N/{N + 2). The leading asymptotics of the spin-spin 
correlation functions at the quantum critical point can be es- 
timated since gi is the spin— 1/2 SU(2)jv primary field with 
scaling dimension 3/2{N + 2).^^ The equal-time spin-spin 
correlation reads then as follows: 



\^a,7i^a^n-\ 



(-l)'"(lnm)i/2 



(14) 



where the logarithmic corrections are the same as for 
spin— 1/2.^^ Similarly, the dimer-dimer correlation function 
is also algebraic with the same power-law behavior but differ- 
ent logarithmic corrections: 



(ea(x)e.(O)) 



(Inx 



1-3/2 



rj,3/(,N + 2} 



(15) 



For N = 2, the semiclassical approach reproduces the 
SU(2)2 transition, first predicted in Ref. 40 using an exact so- 
lution of the perturbation (7) based on the Majorana fermions 
formalism. The location of the transition, obtained within 
our approach, is K/J± = tt^/SG ~ 0.27415. Surprisingly 
enough, this estimate is in good agreement with a previous 
DMRG analysis which found a closure of the spin gap for 
K/Ji_ ~ 0.3 with Jj^ = 1}'^ A second DMRG study locates 
the transition in the window 0.2 < K/J± < 0.26.^^ We will 
come back to this question in our numerical investigation in 
Sec. III. 



2. Jx < 

When J± < 0, i.e. Ai < 0, the minimum condition for 
Vi is now ga+i = ga (a = 1, . . . , A^ — 1). Averaging out 

the g2 N fields, we obtain an effective action on the gi field 

which describes the low-energy properties of the model when 

£:< Ai: 



Seff=NWigi) + Ni-Xi+X2) 



(16) 
When — Ai + A2 > 0, in the semiclassical approach, gi satis- 
fies the constraint Tr(gi) = 0. The effective action (16) be- 
comes equivalent to the non-linear sigma model with a topo- 
logical term 9 = Nn. When N is even, it describes the Hal- 
dane phase of the Heisenberg spin chain S = A^/2 obtained 
when J± — > —00. In contrast, when N is odd, the phase is 
gapless as in the Jj^ > case. 

For sufficiently large negative K, one enters in a new phase 
when — Ai + A2 < which is described by gi — ± J for all A^ . 
This phase corresponds to a uniform dimerization (UD) phase 
where : (Tr(.gi)) = (Tr(.g2)) = . . . = (Tr(.g7v)>- This time, 
the dimerization in a given chain is in-phase with the one in 
flie neighboring chains: ((-l)"Si,„-Si,„+i) = ((-1)"S2,„- 
S2,n+i) = . . . = ((-l)"SAr,„ • S7v,n+i> ^ 0. There is thus 



a quantum phase transition between the Haldane phase (and 
respectively the gapless phase) and the UD phase when N is 
even (respectively odd). The nature of the transition at Ai = 
A2 for model ( 1 ) is the same as in the J± > case and belongs 
to the SU(2)jv universality class. The leading asymptotics of 
the correlation functions at the transition are also given by 
Eqs. (14, 15). 



C. Abelian bosonization approach 

Using a simple semiclassical approach, we have seen that 
an exotic quantum phase transition in the SU(2)7v universality 
class might emerge in model (1) for all N. Here, we present an 
alternative approach, based on Abelian bosonization, which 
confirms the semiclassical prediction. 

In the Abelian bosonization approach, the SU(2)i criticality 
of the decoupled spin— 1/2 Heisenberg chain is described by 
a bosonic field. "*' Introducing N bosonic field $a, the low- 
energy Hamiltonian in the absence of interactions is: 

N „ 

^0 = i E / ^^ ((a.$a)' + (a,$a)') , (17) 

a=l-' 

where e^ {Qa = ^aL - $afl) are dual fields, ^aL.R being the 
chiral components of the Bose fields $a ~ ^aL + ^aR- The 
staggered magnetization (3) and the dimerization operator (5) 
can be expressed in terms of these bosons:"*''^^ 



A 

Trao 

Trao 



■ cos 



V27r Qa) , sin 



2TTQa), 



27r$„ 



(18) 



Using these results, the leading contribution (7) can be 
bosonized: 



n = 



N 
a=l 



dx 



(d^^ar + id^QaY 



dx Y, [91- cos (\/2^ {Qa+l - Qa)) 



g± sill 
34 cos 



27r$a+i ) sm 
27r$a+i ) cos 



27r$a 
2^$„ 



(19) 



where gj^ = JiA'^/aoTr^ and 54 = 36K X^ / gott"^ ■ This model 
has an interesting symmetry content which is explicit within 
the bosonization formalism. Indeed, one can observe that 
model (19) is invariant under the Gaussian duality symmetry 
<1>Q o 8a by fine-tuning the interaction g± = 34: 



■HsD 



N 

is 

a=l 



dx 



{d.^af + {d.QaY 



» JV-1 

■9-L j dx ^ Ices f V27r(6a+1 - ©a) 
•^ a=l 



COS 



(y2^($a+l-$a) 



(20) 



The fine-tuning g± = 54, i.e. K/ J±^ = 7r^/36, corresponds 
to the location of the quantum phase transition for all signs 
of Jl obtained within the semiclassical approach. The low- 
energy effective theory (20), which enjoys a self-dual symme- 
try, has been found in totally different contexts. On the one 
hand, the deconfining phase transition of the 2h-1 -dimensional 
SU( A^) Georgi-Glashow model is controlled by model (20).^-' 
On the other hand, as shown recently, it describes the quantum 
phase transition in dipolar quantum Bose gas."^^'"*' The model 
is also relevant to the competition between superconductivity 
and charge-densy wave or between superfluidity and solid, or 
supersolids in quasi-dimensional systems.^"*"* 

Though the interaction in model (20) is strongly relevant 
with scaling dimension A = 1, a quantum critical behavior, 
stemming from the special form of the interaction, is expected. 
Indeed, model (20) is invariant under two independent global 
U(l) transformations: 



<^a^^a+a 



(21) 



a,/3 being real numbers. The transformation (21) gives a 
U(1)l X U(l)ji' global continuous symmetry of model (20). 
One thus expects the model to be gapless with at least one 
bosonic field protected by the symmetry (21). In fact, it has 
been shown that model (20) displays an extended SU(2)Ar 
quantum critical behavior "^^'^"^ One way to establish this re- 
sult is to use the following conformal embedding 

SU(2)i X SU(2)i X • • • X SU(2)i -> SU(2)jv x Q^, (22) 

where Q^ is a discrete conformal field theory (CFT) with cen- 
tral charge cg^ = N - ■iN/{N + 2) = N{N - l)/(iV + 2). 
It has been shown in Refs. 46 and 63 that the self-dual pertur- 
bation of model (20) corresponds to a special primary field of 
the Qn CFT. A spectral gap opens for the discrete Q^ degrees 
of freedom, leaving the SU(2)jv ones intact. We thus conclude 
that model (20) displays critical properties in the SU(2)Ar uni- 
versality class. The sign of the coupling constant g±^ does not 
play a crucial role for the emergence of this quantum critical- 
ity since one can freely change the sign of the perturbation 
(20) by the transformation: 



$„ 



62 



$„ 



e. 



$ 



2ri+l 



$„ 



e 



2n+l 



e 



2n+l- 



(23) 



We thus confirm the conclusion of the semiclassical ana- 
lysis: interaction (7) exhibits an exotic quantum phase tran- 
sition in the SU(2)Ar universality class for all signs of J^. 
The connection between the two approaches stems from the 
bosonization of the WZNW field ga, which follows from the 
identification (18): 



— iv^STT^a 



9a 




_iV27r$„ 



(24) 



The self-dual sine Gordon perturbation of Eq. (20) is then eas- 
ily shown to be equal to Vi ( 1 2) which governs the quantum 



phase transition in the semiclassical approach. The main in- 
terest of the approach based on the conformal embedding (22) 
is to show non-perturbatively that Vi gives a mass gap for 
the discrete degrees of freedom that is independent from the 
SU(2)Ar ones. 

Though we have shown that the critical properties of model 
(7) are governed by an SU(2)Ar CFT, we are not guaranteed 
that the quantum phase transition of the initial model (1) be- 
longs to the SU(2)jv universality class. Indeed, in general, 
such fixed points with extended criticality are fragile since 
they can be destabilized by several relevant primary operators. 
In addition, the marginal current-current interaction (8), that 
we have neglected so far, might destroy the criticality of the 
transition into a first-order phase transition. The spectrum of 
the SU(2)jv CFT is well-known.^'' The primary operators of 
SU(2)Ar CFT are labeled by their spins j = 1/2, 1, . . . , N/2 
and have scaling dimension 2j{j + 1)/{N + 2). As the con- 
sequence of the SU(2) symmetry of the lattice model (1), 
they take the form (Tr(G))™, m = 2j, G being the SU(2)Ar 
WZNW field. Furthermore, the translational invariance re- 
stricts m to be even since G ^ —G under Tag- The most 
relevant operator is thus the spin— 1 operator, i.e. (Tr(G))^, 
which, according to the semiclassical approach, leads to the 
different phases of the problem. This operator is killed by the 
fine-tuning of the four-spin exchange interaction. The next 
relevant operator is the spin— 2 primary field (m ~ 4) with 
scaling dimension 12/(A^+2) which is a strongly relevant per- 
turbation when A^ > 4. This operator is not forbidden by any 
symmetry of the original lattice. Though it requires a proof, it 
is thus likely that the phase transition is first-order for A^ > 4. 
In this respect, the quantum Monte-Carlo approach of Ref. 
47 points out a first-order phase transition in the N identical 
tubes model with polar molecules when TV > 4. For iV < 4, 
the only possible source of a mass gap is a marginal pertur- 
bation. In this respect, since the SU(2)7v fixed point has been 
obtained for all sign of J±, this criticality is then expected to 
be stable under the current-current interaction at least when 
J_L < 0. We thus conclude that the quantum phase transition 
in the SU(2)Ar universality class occurs for N = 2,3,4 with 
respectively central charge c = 3/2,9/5,2. Large-scale nu- 
merical simulations are clearly call for to reveal the emergence 
of this exotic quantum-critical behavior. 



III. NUMERICAL SIMULATIONS 

In this section, we present extensive numerical simulations 
of the A^-leg spin— 1/2 ladder with cyclic exchange model (1). 
We take the antiferromagnetic leg coupling J|| = 1 as the unit 
of energy. Most data are obtained with DMRG simulations 
with open boundary conditions (OBC). Note that, since a SD 
phase can occur when Jl > 0, we have added two extra sites 
to some of the chains in order to prevent a ground-state de- 
generacy.^^ Typically, we have used 16 sweeps and kept up to 
1600 states which is sufficient for convergence, i.e. having a 
discarded weight of at most 10^®. When K = 0, the physics 
of these spin— 1/2 ladders is well-known*'^ and depends on 
the parity of A^: for N odd (respectively even) the system is 



gapless (respectively gapped), for both signs of J±. When the 
ratio K/J± increases, our numerical simulations confirm the 
existence of a quantum phase transition towards a dimerized 
state that breaks translation symmetry and corresponds to a 
uniform (respectively staggered) pattern between chains for 
J± and K negative (respectively positive). 

Let us mention that Exact Diagonalization (ED) using 
Lanczos technique can also be used to locate the phase tran- 
sition, as there is a clear change of quantum numbers of the 
first excited state across it. This so-called level spectroscopy 
technique has often been used in similar context. ^^•'''^"^*' There 
is also the possibility to investigate the level crossing using 
twisted boundary conditions (TBC), as was done for instance 
in Ref. 71 with a similar model. It turns out that, when avail- 
able, the latter technique tends to reduce the finite-size ef- 
fects in the location of the critical point. Indeed, while peri- 
odic boundary conditions (PBC) engender discrete momenta, 
which can spoil the measure of momentum-dependent opera- 
tors, adding a twist in the boundary conditions allows a con- 
tinuous definition of the momenta and overcomes this prob- 
lem.^- 

Next, we present our numerical results obtained for the 
values of N for which the low-energy study Sec. II pre- 
dicts a phase transition in the SU(2)jv universality class when 
A^ = 2,3,4. Both ferromagnetic and antiferromagnetic J± 
couplings will be considered. 



A. Two-leg ladders (iV = 2) 

1. N = 2 antiferromagnetic rung coupling case: J± — 1.0 

First, we use ED simulations with TBC in order to locate 
the critical point. Using quantum numbers to label respec- 
tively (i) the reflection perpendicular to the rungs (ii) the re- 
flection perpendicular to the legs (i.e. that cuts the bonds 
across the boundary conditions) and (iii) the spin reversal 
symmetry, a simple analysis shows that RS and SD phases 
are identified in sectors (+1, +1, +1) (+1, —1, —1) respec- 
tively.^"* The results for this level crossing are shown for dif- 
ferent system sizes in Fig. 2, which leads to an estimate of 
the critical value at K/J± ~ 0.26. Note that for compari- 
son, we also plot the level crossing found by comparing the 
lowest triplet and singlet states at momentum tt using PBC. 
Using this criterion, we find a critical value K/J± = 0.33. 
However, since it has stronger finite-size effects, we expect 
this result to be less accurate. Our finding is slightly above 
the estimated K = 0.19 found numerically in Refs. 29 and 
30, but in better agreement with K — 0.23 ± 0.03 found in 
Ref. 26. Surprisingly enough, the estimate found in the weak- 
coupling approach {K ~ 0.27) is very close to the numerical 
result obtained for moderate coupling. 

In principle, ED simulations can also give insight into the 
nature of this phase transition.-^ Indeed, the central charge c 
can be extracted from the finite-size corrections of the ground- 
state energy per rung eo{L) = Eq/L'?'^''^^ 
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FIG. 2. (Color online) N = 2 and Jj_ = 1: crossing point with 
PBC (circles) and TBC (squares) versus system size (see text for 
definitions). The most precise extrapolations are obtained using TBC 
data and lead to a critical point at K = 0.26 (dashed line). 



where v is the velocity, which can be extrapolated from^*: 



v{L)^^(E{q^^)-Eiq = 0) 



(26) 



Numerical data, computed at the critical point K = 0.26, are 
shown in Fig. 3 (a,b) and confirm the 1/i^ corrections to the 
ground-state energy per site. We extracted the value of the 
central charge c = 1.52, which is in excellent agreement with 
the expected value c = 3/2 of the SU(2)2 WZNW model. 

However, as detailed below in the DMRG part, while we are 
measuring the correlation exponent, logarithmic corrections, 
coming from marginal operators, seem to be present and spoil 
the accuracy of some of the computations. In such case, it can 
be simpler to compute the scaling dimension x of the lowest 
excitation of the expected SU(2)2 WZNW model by removing 
explicitly these logarithmic corrections^^-^'-^^ 
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[3AEo{S = 1) + AEoiS ^ 0)] (27) 



where AEq{S = 0, 1) are respectively the singlet and triplet 
gaps. Data are shown in Fig. 3 (c) and while finite-size effects 
are non-monotonic, they are relatively weak for large clusters 
and our numbers are compatible with the expected value x = 
3/8. 

We now turn to DMRG simulations to confirm our findings. 
In order to detect the occurrence of a dimerization pattern, 
we measure the local dimerization at the center in the first 
chain, i.e. d{L/2) — (S^ j^/2 ■ ^i.l/2+i)^ ^nd study its scaling 
versus length L. Indeed, according to the low-energy analysis 
of Sec. II, we expect this order parameter to decay as a power 
law at the transition, while it should converge exponentially to 
zero or to a constant below or above the transition. In Fig. 4, 
we plot this order parameter versus system size on a log-log 
plot and it leads to an accurate determination of the critical 
point at K/J± ~ 0.255, which is in excellent agreement with 
our ED result. 

In order to characterize this critical point, we compute the 
entanglement entropy for such parameters. Indeed it is well- 
known^^ that, for a critical point described by a CFT with cen- 
tral charge c, the von Neumann entropy behaves as 
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FIG. 3. (Color online) N = 2 and J± = 1: finite-size scaling of 
various quantities at the critical point K — 0.26: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
SU(2)2 WNZW primary field (see text for details). Data are compa- 
tible with CFT behavior and extrapolations obtain x = 0.375 and the 
central charge c = 1.52. 




FIG. 4. (Color online) N = 2 and J± — 1: dimerization at the center 
f or Af = 2 for various values of K. The critical point coiTesponds to 
a power-law decay (see fit, dashed line) for K = 0.255 (black). 



where d(('\L) = {L/tt) sm{Tri/L) is the conformal distance 
and OBC are used. Fig. 5 shows that our data can be fitted with 
c = 1.50, which is in perfect agreement with the expected 
value c = 3/2 for an SU(2)2 WZNW CFT. 



2. A'^ = 2 ferromagnetic rung coupling case: J± — —1.0 

In the case of a ferromagnetic rung coupling J± = —J\\, 
we need to consider K < in order to induce a transi- 
tion from the Haldane phase to the UD phase. Again, this 
phase transition is detected using ED by comparing the first 
singlet and triplet excited states and plotting this crossing 
value versus system size, as shown in Fig. 6. However, since 
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FIG. 5. (Color online) TV = 2 and Jx = 1: von Neumann 
entropy versus conformal distance d(x\V) at the critical point for 
K = 0.255. 



finite-size effects are quite strong with PBC, we prefer again 
to use TBC and compare eigenstates with quantum numbers 
(+1,+1,-1) (+1,-1, +1) (see definitions in Sec. IIIAl). 
Results for this level crossing are shown in Fig. 6 and TBC 
allow a more precise determination for the critical value at 
K = -0.28. 
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FIG. 6. (Color online) N = 2 and J±_ = —1: crossing point with 
PBC (circles) and TBC (squares) versus system size (see text for 
definitions). The most precise extrapolations are obtained using TBC 
data and lead to a critical point at K = —0.28 (dashed line). 



Similarly to what we have done in the previous section for 
J±^ > 0, we fix this critical value and compute the ground- 
state energies, velocities and gaps in order to determine the 
CFT quantities. Data shown in Fig. 7 allow to determine the 
central charge c = 1.55 and the primary field dimension x = 
0.376 which are in excellent agreement with our prediction of 
a SU(2)2 WZNW universaUty class. 

We now turn to DMRG simulations in order to confirm 
our findings on larger scales. By computing the dimeriza- 
tion at the center (see Fig. 8), we can locate the transition for 
Jf/ J|| ~ —0.275, which is in excellent agreement with our 
ED estimate. At the critical point the dimerization decays as 
a power-law L^" with a = 0.52, which does not agree with 
the expected 3/8 exponent (see Eq. (15)). We believe that 
this disagreement is due to the existence of logarithmic cor- 
rections since a reasonable fit using Eq. (15) can be achieved 
with X ~ 3/8, as also found using ED where we had removed 
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FIG. 7. (Color online) N — 2 and J± = —1: finite-size scaling of 
various quantities at the critical point K = —0.28: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
SU(2)2 WNZW primary field (see text for details). Data are compat- 
ible with CFT behavior and extrapolations obtain x = 0.376 and the 
central charge c = 1.55. 



these logarithmic corrections. 




FIG. 8. (Color online) N — 2 and J± — —1: dimerization at the 
center for various K. Critical point is found at A'/Jy — —0.275 
with a power-law behavior, possibly with logarithmic corrections 
(see text). 

Next we would like to characterize further this critical point 
using the scaling of the von Neumann entropy. However, 
while the use of OBC improves convergence for DMRG, it can 
give rise to so-called Friedel oscillations (see Fig. 9) which 
complicate the analysis. Nevertheless, these oscillations are 
known to originate from the bond modulations,^''**" so we pro- 
pose to use the following fitting form:^' 



SyNii) =A+- lndii\L) + B{S, ■ S,+i) (29) 
6 

where (S; • Sj+i) = J2a(^a,i ■ ^a,i+i) is the sum of the bond 
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contributions. At the critical point, we have found B — —las 
the best parameter, as it completely removes the oscillations 
(see Fig. 10), allowing for a precise determination of c = 1.53 
in excellent agreement with the c = 3/2 analytic prediction. 
Note that we prefer this fitting procedure instead of using PBC 
which also remove oscillations but require much larger nu- 
merical effort. Indeed the number of kept states should be 
much larger for a similar accuracy, typically m^ where m is 
the number of kept states with OBC. Let us also mention that 
other proposals have been made to avoid fitting the oscilla- 
tions,**^ but they only use part of the available data. 
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FIG. 9. (Color online) N — 2 and J± — —1: von Neumann entropy 
versus conformal distance d{x\L) at the critical point K = —0.275. 
Oscillations prevent from a reliable fit of the central charge. 
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FIG. 10. (Color online) N = 2 and J± = —1: same data as in 
Fig. 9 after removing the oscillations due to the bond modulations 
using SvN + {Si ■ Si+i), see Eq. (29). The fit of the central charge 
(dashed line) then gives c = 1.53. 



have checked that the dimerization measured on chain 1 gives 
similar results with the one averaged over all chains. Let us 
also remind that since N is odd, the system is critical in the 
absence of cyclic exchange A' = which implies that the 
transition is now from a gapless phase to a gapped one. 



1. N = 3 antiferromagnetic rung coupling case: J± — 1.0 

In this case we cannot use TBC in ED simulation to get an 
accurate estimate of the critical point. But even with PBC, 
we do not observe any level crossing between the lowest sin- 
glet and triplet with momentum tt (with respect to the ground- 
state) for lengths L up to 12. Therefore we have to determine 
first the critical point using DMRG before performing an ED 
analysis. 

As done before, we measure the dimerization in the mid- 
dle of the first chain. Data are shown in Fig. 1 1 and do indi- 
cate a change of concavity for K/J^^ — 0.37. Note that this 
is not entirely expected since the dimerization should obey a 
power-law in the whole critical phase for K < Kc and con- 
verge exponentially to a finite value in the gapped SD phase 
for K > Kc only. Nevertheless, if we take this heuristic crite- 
rion to detect the quantum phase transition, we get consistent 
results with other approaches (see for instance below for the 
ferromagnetic case where DMRG and ED agree on the loca- 
tion of the critical point). Let us also mention that this change 
of concavity at the transition was also observed recently for 
a similar phase transition in spin— 3/2 chain with three-site 
interaction.'*^ 

At the critical point, we can either fit our data with a power- 
law l/L"-^^, but the expected logarithmic corrections (see 
Eq. (15)) could also give a reasonable fit l/{L^-^\\nL\^^'^) 
(see Fig. 11). 
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FIG. 11. (Color online) A'^ = 3 and J± = 1: dimerization at the 
center for various K. Power-law decay allows to locate the critical 
point for K/J^^ = 0.37. 



B. Three-leg ladders (iV = 3) 

We now turn to the three-leg geometry (iV = 3) using a 
similar analysis. Note that since chains are not equivalent, we 



Now that the critical point is located, we can study the 
scaling of the von Neumann entropy, and as before, we re- 
move oscillations by adding the bond energies and plotting 

SvN{i) + (Si • Si+i) versus the conformal length d{x\L) (see 
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Fig. 12). A simple fit gives a central charge c ~ 1.81 which 
is in excellent agreement with the CFT prediction c ~ 9/5 ~ 

1.8. 
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FIG. 12. (Color online) A'^ = 3 and Jx = 1: von Neumann en- 
tropy versus conformal distance d[x\L) at the critical point for vari- 
ous sizes. Again, oscillations were removed (see text). 
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FIG. 13. (Color online) A*' = 3 and Jx = 1: finite-size scaling of 
various quantities at the critical point K = 0.37: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
primary CFT field (see text for details). 



Assuming now the critical point to be at 7^ = 0.37, we can 
perform a similar ED analysis as for N = 2. Data are shown 
in Fig. 13. Unfortunately, strong finite-size effects, as well 
as different scalings required for L = Ap and L = 4p + 2 
ladders, do not allow for accurate extrapolations of cq and c. 
Nevertheless, if we also perform DMRG with PBC on moder- 
ate L (in which case we have kept up to m = 4000 states for 
L up to 28), we can obtain a reasonable fit of eo versus 1/L^ 
(see Fig. 13 (a)). From this, we extracted c = 1.89, which 
is in reasonable agreement with our previous finding and ex- 
pectation. Note that, since we have used variational DMRG 
energies, one could suspect that either increasing the number 
of states kept or extrapolating with the discarded weight, en- 
ergies will be slightly lower for the largest sizes, thus leading 
to a reduced slope and c value. 

However, ED also allows to extract the scaling dimension 
X by getting rid of the logarithmic corrections. Our numerical 
data support the expected x = 3/10 with an excellent accu- 
racy (see Fig. 13 (c)). 



2. A'^ = Z ferromagnetic rung coupling case: J± — —1.0 

When considering ferromagnetic rung coupling, ED simu- 
lations show a level crossing between the lowest singlet and 
triplet excitations at momentum tt with respect to the ground- 
state momentum. This enables a determination of the critical 
point at Kc = —0.284 (see the inset of Fig. 14), which is ex- 
actly the same as previously found for A'^ = 2. By performing 
a finite-size scaling analysis at fixed K — —0.28, one can ex- 
tract CFT quantities as explained above. Fig. 14 (a-c) present 
finite-size scaling behavior of eg, v and x. As in the previous 
case, size limitation does not allow an accurate estimate of the 
slope of the energy, but still we note that the extrapolation of 
x{L) is rather smooth and leads to x = 0.3 as predicted. In 



order to get an estimate of the central charge, we have also 
performed DMRG with PBC (up to L = 24) in order to com- 
pute eg. Fitting eo and v allows to extract an estimate c = 1.90 
which is not far from the expected 9/5. 
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FIG. 14. (Color online) N = 3 and J± = —1: finite-size scaling of 
various quantities at the critical point K = —0.28: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
primary CFT field (see text for details). Inset: scaling of the level 
crossing versus system size. 

Turning now to DMRG with OBC, the quantum phase tran- 
sition from the critical phase to the UD phase can be located 
as usual on Fig. 15. From the power-law decay of the dimer- 
ization, we determined the transition at K/Jn = —0.275, a 
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value that is extremely close from our ED estimate, thus giv- 
ing us confidence in its accuracy. Indeed, we remind that since 
dimerization should obey a power-law in the whole critical 
phase, it is not obvious that there is a change of concavity at 
the transition. Nevertheless, if we try to fit the dimerization 
decay at the transition, we get exactly the same behavior as in 
the antiferromagnetic case, namely that we can fit either with a 
power-law i^"'*^ (exact same exponent as for J±, K > 0) or 
with the expected behavior from Eq. (15) l/(i"'^| InLj^/"^). 



•X3 



0.01 



0.001, 




FIG. 15. (Color online) A^ = 3 and J± = —1: dimerization at the 
center for various K. The critical point corresponds to a power-law 
decay for iC/ Jy = -0.275. 



Since the von Neumann entropy at the transition exhibits 
large oscillations (see details in Sec. Ill A 1), we use the same 
technique as before and add the total bond energies. In Fig. 16, 
we do observe that oscillations have disappeared and a rather 
good fit can be performed. We obtain c = 1.73, in good agree- 
ment with the expected c = 9/5 = 1.8. 
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FIG. 16. (Color online) TV = 3 and Jx = —1: von Neumann en- 
tropy plus bond energy versus conformal distance d{x\L') at the crit- 
ical point K = -0.275 for L = 128. 



C. iV = 4 

1. N — 4 antiferromagnetic rung coupling case: J± = 1.0 

We start by performing ED simulations and using the cross- 
ing point either with PBC or TBC similarly to the N ~ 2 case 
(see Sec. Ill A 1 for details). Results are shown in Fig. 17 and 
do not allow for a precise determination of Kc which we esti- 
mate to be between 0.3 and 0.5. Therefore we cannot perform 
a CFT analysis at this stage. 
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FIG. 17. (Color online) A'^ = 4 and J± — 1: Crossing point with 
PBC and TBC versus system size (see text for definitions). Ex- 
trapolations are difficult but point towards a critical value between 
0.3 < Kc < 0.5. 

By performing DMRG simulations and plotting the dimer- 
ization at the center (see Fig. 18), we can locate the transition 
precisely for A'/ J|| ~ 0.4. Note this case is similar to iV = 2 
in the sense that both phases are gapped and a change of con- 
cavity is expected at the transition. 

At the critical point the dimerization decays as a power-law 
L~" with a ~ 0.5, which does not agree with the expected 
3/8 exponent. We believe that this disagreement is due to the 
existence of logarithmic corrections as had been observed for 
other values of N. 
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FIG. 18. (Color online) A'^ = 4 and Jx = 1: dimerization at the cen- 
ter for various K. Critical point corresponds to a power-law decay 
for _ft:/ J|| ~ 0.4 



At this critical point, by keeping a large number of states 
(m = 4800 for L — 64), we are able to get a very nice entropy 
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plot (see Fig. 19) and extract c 
with low-energy prediction. 



2.0 in perfect agreement 



2. A'^ = 4 ferromagnetic rung coupling case: J± = —1.0 
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FIG. 19. (Color online) A'^ = 4 and Jx = 1: von Neumann entropy 
versus conformal distance d{x\V) at the critical point K — 0.4 for 
various sizes. 



Going back now to ED simulations at K ~ 0.4 (up to 4 x 10 
ladder), we can perform a similar analysis as before. Data are 
shown in Fig. 20. Unfortunately, strong finite-size effects in- 
hibit us for accurate extrapolations of eg, ?j or c. Similarly, the 
scaling dimension x does not show a sufficiently smooth be- 
haviour to extract an accurate estimate, but it is reasonable to 
extrapolate it between 0.2 and 0.3, which could be compatible 
with the expected 1/4 value. Clearly, we are severely limited 
in system size and it would be useful to be able to remove 
logarithmic corrections in DMRG simulations. 
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FIG. 20. (Color online) A'^ = 4 and Jx = 1: finite-size scaling 
of various quantities at the critical point K — 0.4: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
primary CFT field (see text for details). 




FIG. 21. (Color online) TV = 4 and J± — —1: crossing point with 
PBC and TBC versus system size (see text for definitions). Both PBC 
and TBC allow for an accurate extrapolation to Kc ~ —0.28. 

ED crossing results are shown in Fig. 21 and the use of 
either PBC or TBC data allows to determine the critical point 
K ~ —0.28, again identical to what we had found for N ~ 2 
and A^ = 3. By fixing this value and analyzing our spectra 
(using up to 4 X 10 ladders), we can perform a CFT analysis 
similarly to what has been done for the other N. Data are 
shown in Fig. 22. As for iV = 3, size limitation does not 
allow for a determination of c, however the behavior of x{L) 
is smoother (especially in contrast to the case J± = 1) and 
leads to an estimate x = 0.246 in excellent agreement with 
the expected value 1/4 for /c = 4 SU(2) WZNW. 
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FIG. 22. (Color online) A" = 4 and J± = —1: finite-size scaling of 
various quantities at the critical point K = —0.28: (a) ground-state 
energy per rung eo, (b) velocity v, (c) scaling dimension x of the 
primary CFT field. 

As usual, we can also compute the dimerization at the cen- 
ter by DMRG (see Fig. 23), which allows us to locate the 
transition precisely at K = —0.272, in excellent agreement 
with our ED finding. In particular, it confirms that this critical 
value seems independent from A^ . Note that, at the transition, 
a simple power-law fit would lead to L^"^^ which is not in 
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agreement with the expected x = 0.25; however, taking into 
account logarithmic corrections, we can perform a quite satis- 
factory L-0-25/| InL^/* fit. 
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FIG. 23. (Color online) A^ = 4 and J± = —1: dimerization at the 
center for various K. Critical point corresponds to K/J^^ — —0.272 
with a power-law decay with exponent x = 1/4 up to logarithmic 
corrections. 

Playing again with fitting the von Neumann entropy by re- 
moving oscillations (simply by adding the total bond energies, 
see Fig. 24) allows to get an excellent value for c in agreement 
with c = 2.0. 




d(x\L) 

FIG. 24. (Color online) A^ = 4 and J± = —1: von Neumann en- 
tropy plus bond energy versus conformal distance d{x\L) at the crit- 
ical point K = —0.272 for various sizes. 



IV. CONCLUDING REMARKS 

In summary, we have investigated the nature of the quantum 
phase transitions of A^-leg spin ladders with a cyclic four-spin 
exchange interaction by means of complementary techniques. 
SD and UD dimerized phases emerge from this ring exchange 
on top of the conventional phases of the iV-leg spin ladders. 



Several quantum phase transitions have been studied here de- 
pending on the sign of the ring-exchange interaction K and 
the parity of A^. When TV is even, we can investigate the 
transition between the RS (respectively Haldane) phase and 
the SD (respectively UD) phase when A' > (respectively 
K < 0). In contrast, when N is odd, a transition occurs be- 
tween the c = 1 Heisenberg gapless phase and the dimerized 
phases (SD and UD depending on the sign of K). 

Using a low-energy approach, when the chains are weakly 
coupled, we find that all these quantum phase transitions 
are described by an SU(2)Ar CFT with central charge c = 
3N/{N+2). This result generalizes the findings of Ref. 40 for 
N — 2 where a SU(2)2 emerging quantum criticality has been 
revealed by means of a Majorana-fermion approach, which is 
very specific to the N ^ 2 case. In this respect, we have 
shown that the quantum phase transitions, driven by four-spin 
exchange interactions, belong to the SU(2)Ar universality class 
when iV = 2, 3, 4. For A^ > 4, a relevant perturbation is likely 
to be generated and to drive the system away from criticality. 
A first-order phase transition is then expected for iV > 4 al- 
though further work is clearly called for to confirm or infirm 
this result. 

In close parallel to this field-theory analysis, we have inves- 
tigated numerically the phase transitions for N = 2, 3, 4 by 
means of ED and DMRG approaches. As we have seen, these 
techniques give a precise location of the quantum phase transi- 
tions for all signs of K. The position of the transition does not 
vary too much as a function of N when J± > 0, and seems to 
be independent of N when J± < 0. Surprisingly enough, the 
estimate of the transition obtained within the low-energy ap- 
proach is in rather good agreement with the numerical results. 
The nature of the quantum phase transitions was then obtained 
by extracting numerically the central charge c and the scaling 
dimension x of the lowest primary field of the SU(2)Ar CFT. 
Let us emphasize that while DMRG is the tool of choice to 
measure the central charge c, the presence of logarithmic cor- 
rection prevents a reliable determination of x; on the contrary, 
X can be extracted accurately from ED data as shown in our 
study. In all cases, we found a very good agreement with the 
prediction of the low-energy approach according to which the 
quantum phase transitions should belong to the SU(2)jv uni- 
versality class for iV < 4. In this respect, the A^ = 3 case 
is intriguing due to the prediction of a transition between the 
standard gapless c = 1 phase and the dimerized phases (SD 
and UD phases) that is not of the BKT type with central charge 
c = 1, as for the Ji — J2 Heisenberg spin— 1/2 chain, but 
rather an SU(2)3 transition with the emergence of non-trivial 
critical modes with fractional central charge c = 9/5. 

The field-theory and numerical approaches presented in this 
paper therefore show that four-spin exchange interaction may 
generate an exotic emerging quantum criticality in one dimen- 
sion. In this respect, there is the question whether our analyt- 
ical approach can shed some light on the deconfined quan- 
tum criticality universality class between Neel and UD phases 
in two dimensions.^ Unfortunately, we are not able to deter- 
mine the nature of the phase transition between the RS (or 
the c = 1 gapless phase) and UD phases for A' < when 
N > 2 within our low-energy approach. A phenomenologi- 
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cal approach, describing the 2D system as an infinite array of 
ID spin— 1/2 Heisenberg chains, predicts that the Neel-UD 
transition is governed by an anisotropic 0(4) non-linear sigma 
model with a topological term.**"* Our work is more relevant 
to the investigation of the transition between the Neel phase 
and a staggered VBS which has been observed numerically in 
the spin— 1/2 Heisenberg model with a ring-exchange on the 
square lattice.^"* An estimation of the location of the transition 
by means of ED calculations is K/ J ^ 0.4*^^ which is close 
to our numerical findings in the ladder limit. The nature of the 
universality class of the transition is unknown but it seems to 
be different from the deconfined criticality paradigm.**^ Fur- 
thermore a first-order transition has been found in the Quan- 
tum Monte Carlo simulations of the spin— 1/2 Heisenberg 
model with a special six-spin exchange interaction P Our 
field-theory approach, combined with the ideas presented in 
Ref. 83, seems to suggest that the transition in two dimen- 
sions can be described in terms of an anisotropic 0(4) non- 
linear sigma model. In the future, it will be very meaningful 
to pursue this approach starting from coupled spin chains to 
investigate the nature of two dimensions Neel-SD transition, 
in particular in the easy-plane limit to make connections with 
the phenomenological theory of Ref. 86. 
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Appendix A: Continuum limit 

In this Appendix, we provide the technical details for the 
derivation of the continuum limit of model (1). The ring- 
exchange operator P„+i_„, which performs a cyclic permuta- 
tion of the spins Sa,„, Sa,n+i, Sa+i,„+i, Sa+i,« between two 
consecutive chains, can be written in terms of quadratic and 
biquadratic products of spins as'^'*^^ 

Pn,n+l + Pn.n+l ~ T ^' ^a.n ' Sq+i,,i + Sq,,i+i • Sa+i,„+i 
+ Sa,n • Sq^„+i + Sa+i^„ • Sa+i.„+i + Sa,n ' Sa+l,n+l 
+ Sa+l,n • Sa,ri+1 + 4(S(j^„ • Sa+l,ri) (Sa,?i+1 ' Sa + i,„-|-i) 
+ 4(Sa,ri • Sa,,i+l)(Sa+l,ri ' Sa+l,Ti + l) 
— 4(Sa^„ • Sa+i^„+i)(Sa+i^„ • Sa,n+l)- (Al) 



The spin ladder with a ring-exchange (1) is thus a particular 
case of a general Hamiltonian corresponding to a ladder with 



quadratic and biquadratic spin-spin interactions: 



N 



^gcn — (•-'^ll + Jl) 2_^ /_^ ^aji • Sq^,i+i 
n a— 1 

N-l 

^r / ^ / ^ ^a.n ' ^a+l,n 
n a— 1 

Jd 2_^ 2-^ \^a,n • Sa+l,n+l + Sq^,i-|-1 • Sq-|-1^„) 
n a—1 

N-1 

'''" / y / , [^a.n ' ^a+l,n) i^a,n+l ' ^a+l,n+lj 
n a—1 
N-1 

Ju 2_^ 2-^ (Sa,n ■ Sa,n+l) (Sa+l,n ' Sa+l,n+l) 
n a—1 

N-1 

Jdd2_^ Z_^ (Sa,n ' Sa+l.n+l) (Sa,ii+1 ' Sa+i^„) (A2) 



with Jrr = Ju = —Jdd = 4:K , and J^. = Jj_ + 2K, Jg = 
K, Jd = K for the spin-ladder with a ring-exchange ( Al). 

In the following, we are going to investigate the contin- 
uum limit of the model (A2) in the weak coupling limit when 
\Je,Jr,Jd, Jrr, Jii, Jdd\ ^ J\\, 1-6. in the vicinity of the de- 
coupling point. The continuum limit was done in the two-leg 
case, i.e. A^ = 2, in Ref. 31, but we believe that the identifica- 
tion of the coupling constants of the effective Hamiltonian are 
not correct. To perform such calculations, one needs a free- 
field representation of the fields of Eq. (2) in terms a bosonic 
field ^a with chiral components ^aR,L'-'^^ 

A 



Ua = 1 cos(V27rea),sin(V27rea),-sin(V27r$a 

Ca = COS ( V27r $Q 

TTflo 

1 



J, 



aR,L 



^"^ " 27:a. 



do^-^a 



R,L 



CXp iVStT ^aL , 





where the bosonic fields are normalized according to: 

($aL W *bL (0))=-^ In Z 

47r 



(A3) 



{'^aRiz)<^mm 



sab 



47r 



■ Inz, 



(A4) 



where z = vt + ix (t being the imaginary time) and 
the boson fields satisfy the following commutation relations: 
[^aR{x),^bL{y)] = iS°-^/4:. Within this AbeHan bosoniza- 
tion procedure, the WZNW field ga of Eqs. (3, 5) expresses 

as: 







(A5) 
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This free-boson representation of the SU(2)i CFT is also where we have neglected all chiral interactions, 
useful to derive the different OPEs which will play a crucial 
role in the continuum limit of the general Hamiltonian (A2): 



Jt{z)4{w) 
J%{-z)4{w) 
Jl {z) n^ (w, w) 
J^ (z) n^ {w, w) 
n" (z, -z) 4 {w) 
n" (z, z) J'^ (w) 



6"^ 



Stt^ (z — w) 



,ag-y Jl M 



S'^P 



8TT-^{z-wy 



+ ie 



27r (z — w) 



2TT (z — w) 



—^n" {z,z)n'^ {w,w) 



S"^e{w,w) e^'^'^n'' {w,w) 

Ani (z — w) A'Ki (z — w) 

Airi (z — w) Ani (z — w) 

Airi (z — w) 47ri (z — w) 

Ani (z — w) A'Ki (z — w) 



2ao|z — w\ 



doo|z — w| V 
+ {-z-wf3l{w)) 



2-K' 



\z — w\ 



ao 



(J"^(Ji-Ji,)K«)) 



(j£4 + jfj^)Ku;) 



a^7 



I ,((z-H^ZH 

ao|z- w| 



(A6) 



where z ~ w = vt ~\- iao and z ~ w = vt ~ iao- 

With all these results at hand, one can investigate the lead- 
ing contribution of the continuum limit of the generalized two- 
leg spin ladder (A2). The calculations are very cumbersome 
and our final result reads as follows; 



iLq 



2ttv 



N 



ao 



g I dx J^i-^'aL+JlR) 
{Jr - 2 Jrf) ^ / dx Ua ■ Ua+l 



3ao {Jrr + Jdd + '^Ju) 



N-1 „ 

"^ dx Caea+l 



TT 

N-1 

ao 

a=l 

1 + 4A4 



For the spin ladder (1) with a ring exchange, we finally get 



'l-ring 



+ 



^-^JdxY.{3l,+3l,) 

a.oJ± X! dxua- Ua+l -\ 1 — X] / '^^ *^'*'^°^ 

a=l -^ ^ — 1 "^ 

N 

J2 dx 3aL ■ Jai? [-7 + 2/C (1 

2 (3A^-] 



N 

ao 

a=l 



N-l 



-X' 



■n 

N-l 



ao 2_^ / ^^ ["^aL ■ Ja+li? + -iaR ' Ja+li] 
o=l •' 



J±+AK - AK 



1 + 4A4 



(A8) 



flO X dx 3aL ■ Jafl [-7 + 2 Jf (l - A") 
a=l-' 

■^ [A {Jdd + Jrr + SJu) + Jrr — Jdd — 3 J«] 

N-1 

y ^ / dx [3aL ■ Ja+1_R + Jafl ' Jo+Il] [Jr + '^Jd 



27r2 



{Jdd — Jrr) 



(A7) 
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